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Abstract 

The notion of density of a finite set is introduced. We prove a 
general theorem of set theory which refines the Gibbs, Bose-Einstein, 
and Pareto distributions as well as the Zipf law. 

Suppose that M.^ is a sequence of finite sets tending as n — > oo to an 
infinite set. Suppose that N(A / l^ > ) is the number of elements in the set A4^ n \ 

The set Ai ^ is said to be p( ■ )-measurable if there exists a smooth convex 
function p(-), called a density function, such that the limit 

, p(N(M&>)) 
hm — 

n->oo P[n) 

is finite. This limit is called the p( ■ )- density of the sequence M {n) of sets. 
Let us present a few examples. 

Example 1 Consider the eigenvalues of the k- dimensional oscillator with 
potential 



U(x) = ^2(uixf), x e R k , 

i=i 
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where the u>i are commensurable: 

-A% + U(x)% = Xi%, %(x) e L 2 (R k ). 

Suppose that N\(Xi) is the number of its eigenvalues not exceeding a given 
positive number X. If X —>■ oo, then the limit 

hm — - — - — = k 

a^oo In A 

coincides with the dimension of the oscillator. 

Example 2 Suppose that F is a compact set and Np(e) is the minimal num- 
ber of sets of diameter at most e needed to cover F. Then Np(e) is ln( ■ )- 
measurable and its density coincides with the metric order of the compact 
set F (see J]J/j. 

Consider the set {.Mi} of nonnegative numbers Ai, A2, • • • , A n and the set 
{M2} of integers 1,2,...,N; N = N(n). Suppose that the set {M2} is 
ln( • )@-commensurable and s is its density: 

\nN 
hm = s. 

n — ^00 in n 

Besides, let X(n) be the arithmetic mean of the ensemble of A^: 

A(n) = -^A, 
1 

Suppose that we are given a number E(n). Consider the following cases: 

(1) e < E(n) < X(n)N, e > 0; 

(2) E(n) > X(n)N. 

Consider the set of mappings of {.M2} onto {.Mi}. Two mappings are 
said to be equivalent if their images are identical. Further, we shall only 
consider nonequivalent mappings and denote them by {A^}. 

Suppose that the sum of elements in {.M2} is equal to = Y17=i Ni and 
the bilinear form of the pair of sets {.Mi} and {.M2} satisfies the condition 



Y^NiX, <E(n) 



i=l 



in case (1) and 



i=i 



> E(n) 



in case (2). 

Note that the set {.M3} is lnln( • )-measurable. 

Without loss of generality, we assume that the real numbers Ai, A2, . . . , A n 
are naturally ordered, i.e., < A« < Aj+i, and we split the interval 1,2, ... ,n 
into k intervals (to within 1), where k is independent of n: 

l,2,...,ni, n 1 + l,n 1 + 2...,n 2 , n 2 + 1, n 2 + 2, . . . , n 3 , 

k 

n k - 1 + l,n k - 1 + 2,...,n k , ^ n t = n; 

1=1 

here I — 1, . . . , k is the number of the interval. 

Denote by A/, I — 1, 2, . . . , k, the nonlinear average of Aj over each interval: 



where ip a /3( x ) is the two-parameter family of functions and <& a p is its inverse: 
$0/9(^0/9 (z)) = 1; namely, 

for s > 1; (1) 
for s — 1; 

for < s < 1. 

The parameters ct and /3 are related to jV(ra) and E(n) by the conditions 

n n 

= N(n), Yl X ^{\) = E(n). (2) 



(a) 


Ipocp = 


ae~ px 


(b) 


*Paf3 = 


1 


aeP x — 1 


(c) 


Ipap = 


1 


(3x + In a 



i=l 



1=1 



Consider the subset .A C 7W 3 : 



^={e(X>-^(A/)) <a|, 

^ 1=1 ' J 



(3) 



where 



{VN\n 1/2+e N for iV < n, 
v^ln 1/2+e n for iV ~ n, 
^L\n 1/2+e n for iV > n 



(4) 



is called the resolving power. 

Theorem 1 The following inequality holds: 

M(M 3 \A) < C C 
N{M 3 ) ~ n k N k) 

where k is arbitrary and C is a constant independent of n and N. 
Proof. Obviously, 

J\f{M 3 \ A} = We{jV(n)i<7(n) -£^W^,"<»> 
{N,}^ ^ i=i J 

E( xe { E(E n *-m\i)) 

{Ni} ^ ^ 1=1 H=n ( _! ' 



(5) 



A 



Here the sum is taken over all integers N iy Q(x) is the Heaviside function, 

Q{x) 



1 for x > 
for x < 0, 



and 5k lt k 2 is the Kronecker delta, 



'ki,k 2 



1 for ki = k 2 , 
for ki 7^ /c2- 



Let us use the integral representations 

„~uN pit 



0~N,N> ~ 

©(y) 



2tt 
1 



ie -iN4> e vN> e iN'<) 



— 7T 
OO 



fife 



X — I 



(6) 
(7) 
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We have 



Denote 



df?e"^ = - . (8) 



{Ni} 1=1 
ni -. 



»=nj-i 



and 

T(E,N) =N{M 3 }. 
Since jV{M 3 }(E) < Af{M 3 }(E + e) for e > 0, we have 

POO 

Z(P, N)>/3 dE' T(E', N)e~ fiE ' = T(E, N)e 

J E 



-f3E 



Therefore, 
But, by (6), 



M{M 3 ) < Z((3,N)e l3E . 



-UN /-7T 



Z{0,N) 



2% 



/7T 



da e- iNa C(p,u + ia); 



hence 



N{M 3 \A} 

-vN+(3E fir 



(9) 
(10) 

(11) 

(12) 



< 



2tt 



x e 



exp(-^iV0) (exp{(-/3^^A.) + (i<f> + u ) Nj 



E( E ^--iMa,)) 

1=1 S="i-i ' 



-A 



where /3 and v — — In a are real parameters for which the series is convergent. 



5 



Estimating the right-hand side, carrying the modulus through the integral 
sign and then through the sign of the sum, and integrating over 0, we obtain 

Af{M%\A} <e- uN expPE(n)J2 ex v{-pit, NiXi + uNi } ( 13 ) 

x0 { E(E N i-M\i)) -a}. 

^ 1=1 \=m^i ' ' 

Let us use the following inequality for the hyperbolic cosine cosh (a;) = (e x + 
e~ x )/2: 



Y[cosh(xi) > 2- k e 5 for all x h ^ N > 5 > °- 



(14) 



i=i 



Hence, for all positive c and A, we have the inequality (cf 013]) 







1=1 \=ni-! ' 



k , n t 

-A J> < 2 fe e - cA JJcoshfc J2 Ni-c^piXi 
i=i ^ 



t=ni_ 1 



We obtain 

A/"{A4 3 \ .A} < 2 k e- cA exp{PE{n) - uN) 



(15) 
(16) 



x 



s n s k , ni 

exp I -P NiXi + vN * \ II cosh ( £ ciVj ~ C ^(A, 

{TV;} ^ i=l ^ i=l H=ni_i 



X 



H + c ' ^ exp(-c^ a/ 3(A/)) + - c, /3) exp(c^ a/3 (Ai))) 



!=1 



Let us apply Taylor's formula to <^(z/ + c, /3). Namely, there exists a 7 < 1 
such that 

ln(0(^ + c, /?)) = In $(,/, /3) + c(ln P) + jQn Ci)> + 7C, /?)■ 



Obviously, 



d 

— InO = ipaA^i)- 



Let c = A/D(v,/3), where D(v,/3) = /3). The right-hand side of 

relation (16) is equal to 

Imposing the following constraint on A: 

D {^ + WJY I} )^ {2 -^ D{ ''^ 

where e > 0, and taking into account the fact that D{y,(3) is monotone 
increasing in v, we finally obtain 

Af(M 3 \A) < 2 k e^e- uN C{y ) l3)e-^ 2 l D ^\ 

Next, let us estimate ((u,(3). 

The following lower bound for Z(/3, N) was obtained in j2], relation (95): 

C(z/, P) < y f TW(yJ)Z{P ) N)e v ' N , (17) 
where v' = v'(/3, N) is determined from the condition 

n 

'£ l Zi{v,P) = N. (18) 

i=i 

Suppose that /3 = /3' is determined from the condition 

n 

Y J Ui{v,P) = E{n). (19) 

i=l 

Since Z(/3,N) is determined by the integral (11), its asymptotics given 
by the saddle-point method (the stationary phase due to Laplace) yields a 
unique saddle point for a = 0. 

The square root of the second derivatives with respect to a will appear 
in the denominator. As a result, we obtain 

< Ce Nl/ 'e^' E ^D(p',p')M{M 3 }, 
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where C is a constant. Therefore, we finally obtain 1 



J\f{M 3 \A} ^ k rn(l/ ari-eA'/DV,?) 



AT{M 3 } 



<2 k CD{u',(3')e-^ /n(u ^K (20) 



Further, it is easy to estimate D(y, /3) as a function of N, n: D ~ N for 
N < n, while, for n 3> N, the estimate for D yields the relation D ~ N 2 /n. 
Hence we obtain the estimate for J\f(Ms \ A) / Af {A4 3} , given in the theorem. 

Example 3 For the case in which s > is sufficiently small (and hence, 
Yli=i Ni = N not very large), the Bose-Einstein distribution is of the form 



ni e -Ai/3 



i=ni_ 1 



(21) 



where (3 = l/(kT), T is the temperature, and k is the Boltzmann constant. 

In the case of a Bose gas, for s < 1, we have a distribution of Gibbs type, 
i.e., the ratio of the number of particles on the Ith interval to the number of 
particles on the (I + l)th interval obeys formula (21). 

Example 4 In the case s > 1, we obtain a refinement of the Zipf-Mandelbrot 
law namely, 



However, if s is close to 1, then it is better to use relation (b) in (1), which 
uniformly passes into relation (c) and relation (a). 

Note that if all the A, on the Ith interval are identical and equal to \Q\ 
then Ni/ ni ~ and since n\ ~ N^ s , it follows that, in this case, we 

obtain the Zipf-Mandelbrot formula. 

Example 5 (relation between the sales volume and the prices on the stock 
market) Let us now consider the relation between the prices and the number 
of sold (bought) shares of some particular company on the stock market. 

Since the number n,i of sold shares of that company during the ith day is 
equal to the number of bought shares and \ is the price of the shares at the 

1 A lower bound for N{Mz] was obtained by G. V. KovaP and the author in 0] without 
recourse to the saddle-point method. 
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end of the day, we set, averaging over ni days, the nonlinear average price 
as 

\=n t ' 

•where <\>{x) = l/(x + u), v = const, and is the function inverse to 4>(x). 
Then, by Theorem 1, we have 

n^A^X), (23) 

where A is a constant. 

Thus, the stock market obeys the refined Zipf-Mandelbrot law if all the 
types of transactions are equiprobable (see 0, 

In conclusion, note that although the theorem is stated in terms of set 
theory owing to the fact that we have introduced the notion of equivalent 
mappings, it belongs, most likely, to number theory. Under the same con- 
ditions, considering the set of mappings of the set M.^ onto the set 
without the condition for the equivalence of mappings, i.e., considering all 
mappings, we can obtain a similar theorem that will only be relevant to the 
refinement of the Gibbs distribution. At the same time, such a theorem 
is related to information theory and a generalization of Shannon's entropy. 
Here the estimate has special features, and the corresponding article will be 
published jointly with G. V. Koval'. 
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